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^ Abstract 



We have calculated inclusive direct CP-asymmetries for charmless i? ± -decays. After summing 
large logarithms to all orders the CP-asymmetries in AS = and AS = 1 decays are found as 

a CP (A5 = 0) = (2.0S)%, a CP {AS = I) = (-1.0 ±0.5)%. 

These results are much larger than previous estimates based on a work without summation of 
large logarithms. We further show that the dominant contribution to acp (AS = 0) is propor- 
tional to sin 7 ■ |Kf>/Ki&|- The constraints on the apex (p,rj) of the unitarity triangle obtained 
from these two CP-asymmetries define circles in the (p, rf) -plane. We have likewise analyzed the 
information on the unitarity triangle obtainable from a measurement of the average non-leptonic 
branching ratios Br(AS = 0), Br(AS = 1) and their sum Bt^l(B — > no charm). These 
CP-conserving quantities define circles centered on the p-axis of the (p, rf) -plane. We expect a 
determination of IV^/V^I from Br NL (B — > no charm) to be promising. Our results contain 
some new QCD corrections enhancing Br(AS = 1), which now exceeds Br(AS = 0) by 
roughly a factor of two. 



'Work supported by BMBF under contract no. 06-TM-874. 

2 e-mail:alenz@MPPMU.MPG.DE 

3 e-mail :nierste @ mail. desy. de 

4 e-mail:Gaby.Osterm aier@feynman.t30.physik.tu-muenchen.de 



2 



1 INTRODUCTION 



1. Introduction 



CP-violation is a litmus test for the Standard Model, which parametrizes all CP-violating quan- 
tities by a single parameter, the complex phase in the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix. The related amplitudes are further suppressed due to the smallness of CKM elements and 
loop graphs, so that new physics effects may become detectable. CP-violating observables are 
commonly expressed in terms of the angles a, f3 and 7 of the unitarity triangle. Yet we can 
determine its shape not only from its angles, but also from the length of its sides, which are ob- 
tained from measurements of CP-conserving quantities. This interplay is a special feature of the 
CKM mechanism. In order to overconstrain the unitarity triangle one must find sufficiently many 
theoretically clean observables. While, for example, f3 can be extracted without hadronic uncer- 
tainties from the mixing-induced CP asymmetry in B d — > J/tpK s , the angle 7 is notoriously 
hard to measure in experiments with B d and B^ mesons. 

Direct CP-violation in exclusive 5 ± -decays does not help to determine any of the angles be- 
cause of the unknown strong phases in the decay amplitudes. On the contrary direct inclusive 
CP-asymmetries can be cleanly predicted, because quark-hadron duality allows the reliable cal- 
culation of strong interaction effects within perturbation theory. Such direct inclusive asymme- 
tries have been analyzed in [1-6] and mixing-induced inclusive CP-asymmetries studied in [7] 
are now investigated by the SLD collaboration [8]. Semi-inclusive direct CP-asymmetries have 
been studied in [9]. While inclusive final states are experimentally difficult to identify, inclusive 
branching ratios are huge compared to exclusive ones. As we will see in the following, inclusive 
CP-asymmetries in charmless decays have a promising size, so that it is worthwile to study them 
experimentally. Further they can be obtained from branching ratios only and therefore do not 
require an asymmetric fi-factory. 

In this paper we calculate direct inclusive CP-asymmetries in charmless S ± -decays extending 
our recent calculation of decay rates in [10]. In [10] the corresponding branching ratios have 
been calculated in renormalization group improved perturbation theory including the dominant 
contributions of the next-to-leading order. In the following section we set up our notations and 
summarize previous work on the subject. In sect. 3 we analyze AS = decays. We discuss 
the relation of the CP-asymmetries to the angles of the unitarity triangle and their impact on 
the determination of the improved Wolfenstein parameters p and rj. Here we also investigate the 
constraint on (p, rj) imposed by a measurement of the average branching ratio of B and B decays 
with AS = 0. In sect. 4 we repeat the procedure for AS = 1 decays. Readers mainly interested 
in phenomenology may draw their attention to sect. 5, where we will give numerical predictions 
for the newly calculated quantities. In sect. 5 we also predict the total charmless non-leptonic 
branching ratio of the B meson and exemplify how the unitarity triangle is constructed from the 
branching ratios and CP-asymmetries. Finally we summarize our findings. An appendix contains 
details of our analytical results. 
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2. Preliminaries 

We start our discussion with B decays corresponding to the quark level transition b 
q = u,d, s, c. They are triggered by the \ AB\ = 1, \ AS\ = hamiltonian H: 



H 



G F f 



^§ | E q {CQ- + CQ-) - £ E CjQj j+h.c. 



v; b v qd 



qqd, 



(1) 



Here Qi'2 are the familiar current-current operators, which originate from the tree-level W- 
exchange in b — > ccrf and 6 — > uud. Further V — {3, . . . 6, 8}, and Q3-6 and Qs are the penguin 
operators. More details can be found in [10], where the numerical values for the Wilson coeffi- 
cients Ci are tabulated. For the following we only have to keep in mind that the coefficients C 3 _ 6 
and C 8 accompanying £ t * are much smaller in magnitude than C\ and C 2 . 
Now we express the decay rate for b — > qqd as 



64tt 3 



Re 



if i 2 r 



+ I6l 2 r« + |e c | 2 r cc + u c *r, 



+ tel^u + 6Cr tc • (2) 



The coefficients encode the various contributions of the different operators in (1). For example 
in b — > uud the interference of the tree diagram of in Fig. 1 with the penguin diagram of Q c 2 
with q' = c in Fig. 2 contributes to T uc . The average branching ratio for the decay of B 
some inclusive final state X reads 



± into 



Br 



T(B + -> X) + r (b- -> x) 



2r 



(3) 



tot 



Similarly we define the CP-asymmetries as 

1 



A 



CP 



Br (B + ^X) - Br (b~ -> X 



CP 



5r 



(4) 



Of course the average branching ratio in (3) may also be considered for B d and B d instead of 
B + and B~. We do not consider small spectator effects in this work, so that all given formulae 
for Br likewise apply to the neutral B mesons. We will classify the inclusive final state X by 
its strangeness quantum number S. Hence if also B s mesons are included in the consideration of 
Br, the strangeness of X must be corrected for the non-zero strangeness of the spectator quark. 
Our strategy is to express Br and the CP-asymmetries in (4) in terms of the T^'s. The con- 
straints for the CKM matrix obtained from measurements of Br and a C p are most conveniently 
expressed in terms of the improved Wolfenstein parameters (p, rj) [11]. The calculation of (p, rf) 
from Br and acp involves certain combinations of the T^ 's, for which we will derive compact 
approximate formulae. The exact expressions for the T^ 's can be found in the appendix. 
The three angles of the unitarity triangle are 



arg(-&0 = 



a. 



arg = 7- 



(5) 
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2 PRELIMINARIES 





Figure 1: Tree diagram of Q\ 2 . 



Figure 2: Penguin diagram involving Q\ . 
It contributes the absorptive part necessary 
for a non-zero acp- 



Then Br and Acp are easily found as 



Br = F ■ 



\V cb ? 



- |£u£c| cos7Rer uc - |£ t £ u | cosaRer tu - |£ t £ c | cos(3ReT tc ] (6) 



CP 



[- Ifufd sin7lmr uc - \£ t £ u \ sinalmT^ + \£ t £ c \ sin/3Imr te ] . (7) 



The common factor F reads 



rjexp 

0.1045 



0.715 + 3.0(x c - 0.3) + U(x c - 0.3)' 



1 -0.04 log 



JL 

m b 



(8) 



Here x c = m c /m b and /x = 0(m b ) is the renormalization scale. F is inverse proportional to the 
total decay rate T tot , which we calculate via T tot = T SL /B SL from the measured semileptonic 
branching ratio B SL . The numerical approximation in (8) holds to an accuracy of 1 % in the 
range 0.25 < x c < 0.35 and for variations of the renormalization scale /i in the range m b /2 < 
/j, < 2m b . The exact expression can be found in the appendix. 

From (7) one can nicely verify that one needs two different CKM structures and a non-zero 
absorptive part ImT^ in order to obtain a non- vanishing Ac p. It is known for long that the 
CPT theorem correlates the CP-asymmetries for different subsets of final states X in (4) [2,5]. 
For example A CP (A\C\ = 2, AS = 0) = -A CP (A\C\ = 0, AS = 0), where (A|C| = 2, 
AS = 0) denotes the decay into the inclusive final state with total strangeness zero containing a 
c and a c quark, while A|C| = corresponds to a charmless final state. In the following we will 
focus on charmless final state and omit "A|C| = 0" in our notation. The non-zero contributions 
to Acp(AS = 0) come from the absorptive parts of penguin diagrams (see Fig. 2) involving 
the annihilation process (c,c) — > (q,q), q = u,d,s. We have illustrated the leading 0(a s ) 
contribution to Imr uc , lmT tc and Imr (u in Figs. 3 and 4. The results of all possible operator 
insertions into Figs. 3 and 4 can be expressed in terms of a single function g(m q > /m b , fi/m b ), 
e.g. Imr tc oc lmg(x c , /J,/m b ) (for details see the appendix and [10]). We will need some special 
values: 



0(0,1) = -0.67- 0.93i, g(x c = 0.3, 1) = -0.69- 0.23i, 0(1,1) = 0.28.(9) 



Figure 3: Diagram constituing Imr MC to order a s for AS = decays. The right cut marks the 
final state uud. The left cut denotes the absorptive part of the penguin diagram of Fig. 2. 



The imaginary part of g is /i-independent. Incidentally we will omit the second argument of g. 
Let us now look at the CP-asymmetry related to a specific quark final state, for definiteness we 
consider uud: The contribution from r uc depicted in Fig. 3 involves Q\ and and is therefore 
proportional to Cf , while T tc and T tu involve £ t * and a small penguin coefficient C 3 _ 6 thereby. 
Now Imr uc oc lmg(x c ) = —0.23 for x c = 0.3. The smallness of Im g(0. 3) compared to Im g(0) 
in (9) reflects the fact that Im g (x c ) vanishes for x c > 1/2. Yet Gerard and Hou [2] have made the 
important observation that this kinematic suppression is absent in the higher order contributions 
to r uc , so that the result of Fig. 3 receives a correction of order a s (mb) /n • Im g(0) /lmg(x c ) ~ 
30%. But these unsuppressed terms cancel in the sum Acp (uud) + Aqp (ssd) + Acp (ddd) = 




A CP (AS = 0) = —A C p (A\C\ = 2, AS = 0), because the latter asymmetry vanishes in the 
kinematically forbidden region x c > 1/2 [2,5]. In this work we will only calculate the inclusive 
CP asymmetries for charmless AS = and AS = 1 decays and therefore do not need to 
include terms of order a 2 . This, however, is not true for the separate inclusive CP-asymmetries 
Acp (ssd') and A C p (ddd'^j, d! = d,s, calculated to order C 2 C^ & a s in [4]. In addition the 
0(a^)-contributions to these quantities involve the large results of "double penguin" diagrams 
proportional to C\ corresponding to the square of the diagram in Fig. 2. 

Still there is an important difference between our calculations and those in [2]: We use the ef- 
fective hamiltonian of (1), while Gerard and Hou perform their calculation in the full theory and 
thereby invoke large logarithms, which are summed to all orders in our approach. These large 
logarithms lead to an apparently large contribution of order o? s in [2], which had been found to 
cancel the leading contribution of order a s numerically, so that the authors of [2] have claimed 
the total inclusive asymmetries to be vanishingly small, of order of a few permille. As we will 
see in the following, the correct summation of the large logarithms leads to a different result: The 
inclusive CP-asymmetries acp(AS = 0) and acp(AS = 1) are sizeable, of the order of two and 
one percent, respectively. 

CP-asymmetries with resummed large logarithms have also been calculated in [6], but for the 
case of a light m t ~ 15 GeV. In [6] therefore no penguin operators Q3-6 and Q 8 appear. The 
actual numerical results for m t ~ 170 GeV are substantially different. In [6] also the observation 
has been made that corrections of order a 2 s are small for a C p(AS = 0) and a C p(AS = 1). 
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3 AS = DECAYS 





Figure 4: The diagrams show the contributions to Imr te and lmT tu to order a s . lmT tc corre- 
sponds to q' = c with the right cut denoting the final state with q = u, d or s. Likewise IinT^ 
is obtained for q' = u, but now also the left cut marks a possible charmless final state. Then the 
right cut denotes the absorptive part of a penguin diagram with internal quarks q = u, d, s or 
c. The contributions of this figure are suppressed with respect to those of Fig. 3, because they 
involve a small penguin coefficient C 3 _ 6 . 



3. AS = decays 

We first look at the dominant contributions to Br and a C p'- Keeping only the lowest nonvanishing 
order in a s and neglecting the contributions of the small penguin coefficients one finds 



Br(AS = 0) = F\V t 



ud\ 



K 



ub 



cb 



a C p(AS = 0) = -lmF v 



V, 



cd 



V, 



ad 



V, 



cb 



V,. 



ab 



sin 7. (10) 



Hence from Br one can determine |T4&/Kb|> because F and V u d are well-known. Likewise acp 
measures the product of sin 7 and | V c b/V u b |. The corrections to (10) stemming from the penguin 
coefficients and higher order corrections to Br are reliably calculable and small. The best way to 
exploit (10) and to include these corrections is the use of the improved Wolfenstein parameters 
A, A, p and 77 [11]. Then Br of (6) reads 



Br(AS = 0) = L B (p 



Po) 2 + t 



K 



(ID 
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Here 



L B = FX 2 [T tt + T uu - Re r tu ] (12a) 
2Tg + Re [r MC - Ttc - T tu ] 

Po ~ 2 [r tt + r uu - Re r tu ] (12b) 



K = 



1 



4 [r tt + r uu -R e r tu ] 



-AT uu T tt + 4T uu Re T tc + [Re T uc } 2 - 2Re T uc Re T 



tu 



- 4r cc r uu + 4r cc Re T tu + [Re T tu ] 2 - 2Re T tu Re r tc 
-4r tt r cc + 4r«Re r uc + [Re T tc f - 2Re T tc Re T uc ] . (12c) 

We stress here that our notation of Br(AS = 0) only comprises non-leptonic decays, but not the 
semileptonic decay B — > X u tp£, which is measured in a different way. In addition to the quark 
final states uud, ssd and ddd we have also included the decay b — > dg, which gives a small 
contribution of order 3% to Br(AS = 0), but has a non-negligible impact on K and p - Notice 
that the Wolfenstein parameter A drops out in (12). The corrections to the formulae in (12) are of 
order A 6 and therefore negligible. From (1 1) one sees that the measurement of the CP-conserving 
quantity Br defines a circle in the (p, rj) -plane centered at (p , 0) with radius R B , where 

Rl - *<^>+A-. (13, 

The center (p , 0) and K are independent of the measured Br, they vanish in the limit considered 
in (10). For the constraint from the CP asymmetry we likewise define 

L a = Im [Ttc - Vtu - . (14) 



uc\ 

r<t + T uu — Re T tu 



Then 



A CP (AS = 0) = L a L B % a CP (AS = 0) = L _ _ V _ 2 -. (15) 

Again the corrections to (15) are suppressed with four powers of A and therefore negligible. Now 
(15) reveals that a measurement of acp likewise fixes a circle in the (p, rf) -plane. This new circle 
is centered at (p , r} ) and its radius equals R a with 

= 2 OCP (A g = 0) - + ^ 

Again in the approximation with K = p = adopted in (10) the circle defined by (15) is cen- 
tered exactly on the 77-axis. Its radius equals fj , so that it passes through the origin. In (10) sin 7 

comes with | Kb/Kb |> which is inverse proportional to ^p 2 + rj 2 . The geometrical construction 
of 7 from a C p corresponding to (10) is therefore done by intersecting the circle in (15) with the 
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4 AS = 1 DECAYS 



one centered at (0, 0) stemming from any measurement of \V u b/V c b\. Of course any other infor- 
mation on the apex (p, rf) of the unitarity triangle can be included in the usual way, and ideally 
the hyperbola from [12, 13], the circle from Am B [13] and the new circles in (11) and (15) 
intersect in the same point (p, rj) — or we may find new physics. 

We close this section by giving compact approximate expressions for the quantities in (12) and 
(14), which enter the circles defined by (12b), (13) and (16): 



L B = (0.0362 + 0.151 (x c - 0.3) + 0.58 (x c - 0. 

.(l -0.12 In +0.02 In 2 


3) 2 ) 


uncertainty: 0.80% 


L a = (0.00734 - 0.0905 (x c - 0.3) + 0.220 (x c - 0.3) 2 ) 

.(l -0.22 In +0.08 In 2 JL) 

v rrib TTib J 


uncertainty: 3.6% 


K = (- 0.0133 + 0.0219 (x c - 0.3) + 0.032 (x c - 0.3) 2 ) 

•(1-0.44 In +0.16 In 2 JL) 

v m b m b J 


uncertainty: 3.3% 


p = (- 0.0254 + 0.034 (x c - 0.3) + 0.12 (x c - 0. 

■(1 + 0.03 ln 2 ^ 

V rrib 1 


3) 2 ) 


uncertainty: 1.8% 



(17) 



In the last column we have listed the error of our approximate formulae compared to the exact 
expressions for the range 0.25 < x c < 0.35 and 0.5 < fi/rrib < 2.0. Further a s [Mz) = 0.118 and 
L B is calculated with B S l = 0.1045. The p-dependence in (17) results from the truncation of the 
perturbation series and is small in L B , for which the dominant next- to-leading order corrections 
are known. A future calculation of the full 0(a s ) corrections to Br and the 0(a 2 s ) corrections to 
Aqp will change the numbers in the first brackets in (17) by a term of order a s (rrib)/n and will 
reduce the size of the coefficients of ln(p/mfe) and ln 2 (p/m b ). 

4. AS = 1 decays 

To obtain the \ AS\ = 1 hamiltonian from (1) we must simply replace £ q by £^ = V* b V qs . Instead 
of (5) we invoke the CKM angles 

arg (-d s) d s) *) = -0 = -( 7 -A 2 ^)+0(A 4 ), 
arg(-e } d s) *) = -e = -X 2 v(l + X\l-p))+0(\ e ), 
arg (-d s) 4t } *) = 7-vr + 0(A 4 ). 

Hence the corresponding unitarity triangle with angles n — 7, and e is squashed. In the limit of 
vanishing penguin coefficients one has 

A CP (AS = 1) oc \V ub V cb \sm-f. 
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Yet an approximate formula for acp (AS = 1) similar to (10) cannot be found, because the 
tree-level contribution to Br (AS = 1) is CKM suppressed and the different T^'s are equally 
important. An analogue of (10) would involve more than one CKM angle. 
Next we express Br, A C p and a C p as in (1 1) and (15): 



B7(AS = 1) = L' B [(p-p> o y+f-K> 
A CP (AS=1) = L' a L' B rj, a CP (AS = 1) = L'„ 



V 



a (p-p' y + v 2 -K>- 



(18) 
(19) 



The primed coefficients read 



L' = 



L' 



K' 



FA 4 (l + A 2 ) (Y uu + Y tt -ReY tu ), 
-2T tt + Re (T tu + T tc - T uc ) 

2A 2 (i + A 2 )(r UM + r tt -R e r tu )' 
im (v uc + v tu - r fc ) 

A 2 (l + A 2 ) (Y uu + Y tt -ReY tu y 
-(l-\ 2 )(T cc + T tt -ReT tc ) 



[-2T tt + Re (T tu + T tc - T uc )f 



A 4 (1 + A 2 ) (T uu + T tt - Re r tu ) 4 A 4 (1 + A 2 ) (T uu + T tt - Re T 



(20) 



tv J 



In L' B , ~p' and K' we have kept corrections of order A 2 and omitted corrections of order A 4 and 
higher in accordance with the adopted improved Wolfenstein approximation [11]. The powers 
of A in the denominators of p' and K' are partially numerically compensated by the smallness 
of the penguin coefficients entering the Tj/s in the numerators. The corresponding approximate 
formulae read 



L' B = (0.00185 + 0.0077 (x c - 0.3) + 0.03 (x c -0.3) 2 ) 
.(l _ .12 In +0.02 In 2 


uncertainty: 1.0% 


L' a = (- 0.144+ 1.77 (x c - 0.3) - 4 (x c - 0.3) 2 ^ 
.fl _ .22 In + 0.08 In 2 JL) 




uncertainty: 4.3% 


K' = (- 5.08 + 8.4 (x c - 0.3) + 12 (x c - 0.3) ; 

•fl-0.441n^- +0.16 In 2 JL) 




uncertainty: 2.9% 


% = (0.498 - 0.67 (x c - 0.3) - 2.4 (x c - 0.3) 

.(1 + 0.03 ln 2 ^) 


2 ) 


uncertainty: 1.8% 



(21) 



Here we emphasize that in (21) we have not only included the final states with quark contents 
uus, dds and sss, but also the decay b — > s g, which gives a non-negligible contribution to 
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5 PHENOMENOLOGY 



Br(AS = 1) in (18). Further we had to include the contributions to the decay rate stemming 
from the square of the penguin diagram in Fig. 2. These contributions are of order a 2 s , but are 
proportional to Cf and the fourth power of A. They belong to T cc in (2) and amount to 13 % 
of Br (AS = 1). The large contributions of penguin operators and penguin diagrams imply 
that Br(AS = 1) is quite insensitive to p and rj. This is reflected by the large value of K' in 

(21) . Consequently Br(AS = 1) becomes only a useful observable to constrain (p, rj) once its 
experimental accuracy is better than 10 %. 

The geometrical constructions of the circles obtained from Br(AS = 1) and a C p(AS = 1) is 
done in a completely analogous way to sect. 3. One merely has to replace the unprimed quantities 
in (13) and (16) by the primed ones of (21) to obtain the AS = 1 parameters R' B , rj' and R' a . 
Since the denominator of a C p(AS = 1) in (15) depends very weakly on p and rj, a C p{AS = 1) 
is almost proportional to rj and both radius R' a and offset rj' of the corresponding circle are very 
large. This is very different from the situation in AS = decays. 
Finally we mention that 

A CP (AS = 1) = -A CP (AS = 0) (22) 

for m s = 0. This is a consequence of the CKM mechanism of CP violation. The relation in (22) 
receives corrections by terms of order m 2 /m 2 and m s /m b ■ a s (m b )/ir. A larger deviation from 

(22) would be an experimental sign of non-standard CP- violation outside the quark mass matrix. 

5. Phenomenology 

In this section we give numerical predictions for the branching ratios and CP asymmetries and ex- 
emplify, how the apex (p, rj) is constructed from future measurements of Br NL (B — > no charm), 
a C p(AS = 0) anda CP (AS = 1). 

First we express Br NL (B — > no charm) analogously to (11) and (18): 

W NL (B -> no charm) = Br (AS = 0) + Br (AS = 1) = L B (p 2 + rj 2 -K) . (23) 

There is no dependence on 7 here, i.e. p = 0, for the same reason as (22). It is easy to relate L B 
and K to L B , L' B , K and K' . The approximate formulae read 



L B = (0.0380 + 0.158 (x c - 0.3) + 0.6 (x c -0.3) 2 ) 
• fl-0.12 ln^- +0.02 In 2 JL) 


uncertainty: 0.80% 


K = (- 0.272 + 0.46 (x c - 0.3) + 0.7 (x c - 0.3) 2 ) 
•(1-0.42 ln^- +0.15 In 2 JL) 


uncertainty: 2.7%. 



(24) 



As usual the last column lists the error of the approximate formulae for the range 0.25 < x c < 
0.35 and 0.5 < p/m b < 2.0 with a s (M z ) = 0.118 and B SL = 0.1045. 



5 . 1 Numerical predicitions 
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V ^ b =0.06 

Vcb 


^ =0.07 

Vcb 


^ =0.08 

Vcb 


^ =0.09 


=0.10 

Vcb 


0127 +0 - 0057 

U.UIZ/ _Q.0034 


nmu +0.0058 

0.0136 „ .oo35 


0147 +a0060 

U.Ulf/ _o.0036 


015Q +0 0062 

U.Ul^y _o.0038 


1 72 +0 - 0064 

U.Ul 1 L _o.0040 



Table 1: The total nonleptonic charmless branching ratio Br NL (B — > no charm) as a function of 
It is independent of 7. 



5.1. Numerical predicitions 

Next we predict the average branching ratios and the CP asymmetries as a function of iKb/K&l 
and 7. For this we recall the relation of these quantities to the improved Wolfenstein parameters 
(p,77) [11]: 



V. 



ub 



cb 




A Jp 2 + r/ 2 , tan 7 = \ (25) 

v P 



The predictions for the branching ratios can be found in Tabs. 1, 2 and 3. Then acp is tabulated 
in Tabs. 4 and 5. The range of 7 in the tables is the one favoured by the standard next-to-leading 
order [12, 13] analysis of the unitarity triangle from e K and Am B . The central values in the tables 
correspond to the following set of input parameters: 

x c = 0.29, p = m b = 4.8 GeV, 

a s (M z ) = 0.118, m t (m t ) = 168 GeV, B e s f = 0.1045. (26) 

Here m b is the one-loop pole mass. The errors in the tables correspond to a variation of x c = 
m c /m b and the renormalization scale /x within the range 

0.25 < x c < 0.33, 0.5 < fi/m h < 2.0. 

The corresponding error bars are added in quadrature. The experimental uncertainty in a s has a 
smaller impact on the listed quantities, the errors of the remaining input quantities in (26) have a 
negligible influence. 

From a comparison of Tab. 3 with Tab. 2 one realizes that charmless non-leptonic fi-decays 
occur preferably with AS = 1, with Br (AS = 1) exceeding Br (AS = 0) by roughly a factor 
of two: 

- 0.50 ±0.12 for §^ = 0.08. 

Br(A5' = l) \V cb \ 

Most of the dependence on x c stems from the normalization factor F and cancels in ratios of 
different Br's. The yU-dependence of Br (AS = 1) is much larger than the one of Br (AS = 0) 
leading to larger error bars in Tab. 2. This comes from the penguin dominance of Br (AS = 1) 
and the fact that current-current type radiative corrections to penguin operators have not been 
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5 PHENOMENOLOGY 
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Table 2: The average inclusive branching ratio into nonleptonic final states with zero strangeness, 
Br (AS = 0), as a function of |Kb/Kfe| and 7. 



calculated yet. The newly calculated contributions enhance Br (AS = 1) explaining the increase 
of Br NL (B no charm) in Tab. 1 compared to the result in [10]. In order to obtain the total 
charmless branching ratio Br(B — > no charm) one must add twice the charmless semileptonic 
branching ratio Br(B — > X u £v e ), for £ = e and t — /i [14]: 

Br(B^XJV e ) = (o.0012t°££) • (^F^ ■ 
Hence for the input of (26) one finds from Tab. 1 : 



2 



Wr{B - no ctann) = 0.00971*™ + (0.0073^) 

The present experimental result for the total charmless branching ratio reads 

~Br exp (B — > no charm) = 0.002 ± 0.041, 
obtained in [15] from CLEO data [16]. 

We conclude that the measurement of Br (AS = 0) provides a competitive method to deter- 
mine IK^/Vcbl compared to the standard analysis from semileptonic decays. Once a complete 
next-to-leading order calculation is done for the AS = 1 decays, the error bars in Tab. 1 will 
reduce significantly and Br NL (B — > no charm) will likewise become a promising observable to 
measure \ V uh /V ch \. 

The most important results of our calculations, however, are those listed in Tab. 4 and Tab. 5. 
Adding the errors stemming from the uncertainties in IV^/V^I and 7 in quadrature to the ones 
already included in the tables, we predict: 

a CP (AS = 0) = (2.0^;o)%, a CP (AS = l) = (-1.0 ±0.5)%. (27) 



5.2 Construction of (p, 77) 
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Table 3: The average inclusive branching ratio into nonleptonic final states with strangeness one, 
Br (AS = 1), as a function of (Kb/Kfel and 7. 



These results have to be contrasted with those of Table 1 in [2], where predictions for the a C p's 
are given, which are five times smaller than those in (27). This discrepancy is partly due to the fact 
that we sum large logs to all orders whereas this has not been done in [2]. It is further related to 
the use of an extremely small | sin 7 • V c b/V u b\ in [2]. The reduction of the ^-dependence in Tab. 4 
and Tab. 5 requires the calculation of Im T to order a 2 s . The corresponding diagrams are obtained 
by dressing Fig. 3 and Fig. 4 with an extra gluon. A part of this calculation has been performed 
in [3]. In a perfect experiment the detection of a C p (AS = 0) = 2% with ~Br (AS = 0) = 5-10~ 3 
at the 3a level requires the production of 4.5 • 10 6 B ± mesons. This should be worth looking at 
by our experimental colleagues. Finally we remark that our results satisfy 

a C p(AS = 1) • Br (AS = 1) = -a CP (AS = 0) • Br (AS = 0) . 

as required by (22). 

5.2. Construction of (p, rf) 

In this section we exemplify how the circles in the (p, 77) -plane will be constructed from a fu- 
ture measurement of Br NL (B — > no charm), Br (AS = 0), Br (AS = 1), a C p (AS = 0) and 
a C P (AS-\). 

We first show this construction for the CP-conserving quantities. We assume that the three charm- 
less non-leptonic branching ratios are measured as 

Br ml (B — > no charm) = 1.47%, 
W(AS = 0) = 0.50%, W(AS = 1) = 0.97%. 

For illustration we assume an experimental error of 5 % in all quantities and neglect the present 
theoretical uncertainty here by setting // = m b and x c = 0.29. The three circles are defined by 
(11), (18) and (23). To draw the circles we must only read off the coefficients from (17), (21) 
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5 PHENOMENOLOGY 
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Table 4: The inclusive indirect CP-asymmetry for charmless final states with zero strangeness, 
a CP (AS = 0). 



or (24) and calculate the radii Rb, R'b an d Rb from (13). The results are shown in Fig. 5. The 
figure reveals that Br (AS = 0) is a very good observable for the phenomenology of the unitarity 
triangle. This remains true even if the actual theoretical uncertainty of 20 % is included. By 
contrast Br (AS = 1) is not very sensitive to (p, rj) and thereby yields a much poorer information 
on the unitarity triangle. Still the center (p , 0) of the circle largely deviates from the origin, so 
that upper or lower bounds on Br (AS = 1) could help to exclude a part of the (p, ?j)-plane 
allowed by other observables. Also Fig. 5 shows that Btnl{B — * no charm) is a very useful 

observable to determine y^p 2 + if, once the large p-dependence of the entries in Tab. 1 is reduced 
by a complete next-to-leading order calculation of the AS = 1 decay rates. 

The circles from the CP-asymmetries are likewise obtained from (16). Here the measured value 
of a C p enters the ^-coordinate rj or rj' of the center of the circle and its radius R a or R' a . For the 
CP-asymmetries we assume an experimental precision of 20 % and 

a C p (AS = 0) = 2.0%, a C p (AS = 1) = -1.0%. 

The results are displayed in Fig. 6. If one switches off the effects of penguin operators, the circle 
from a C p (AS = 0) touches the p-axis in the point (0, 0). The distance of the points on the circle 

to the origin is therefore proportional to sin 7, so that a C p (AS = 0) measures sin7/y / p 2 + rf 
in this limit as found in (10). The circle from a C p (AS = 1), however, looks totally different: rj 
and R' a are so large that only a small fraction of the circle can be seen in Fig. 6. acp (AS = 1) 
weakly depends on p and yields good information on rj. Hence from Fig. 6 we learn that inclusive 
CP-asymmetries yield interesting information on the unitarity triangle, which is complementary 
to the one obtained from other observables in the B system. Alternatively one can multiply a C p 
with the measured Br and obtain A C p of (4), which defines a horizontal straight line in the 
(p, rp)-plane (see (15) and (19)). 
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Table 5: The inclusive indirect CP-asymmetry for charmless final states with strangeness one, 
a CP (AS = l). 



6. Ten messages from this work 

1) Inclusive direct CP-asymmetries in charmless 5 ± -decays are larger than previously be- 
lieved: 

a CP (AS = 0) = (2.0+1:0)%, a CP (AS = l) = (-1.0 ±0.5)%. 

2) The dominant contribution to a C p (AS = 0) satisfies 

sin 'v 

a CP (AS = 0) oc 7 (28) 

I VW vcfil 

with small and calculable corrections. 

3) The constraints on the apex (p, rf) of the unitarity triangle obtained from a measurement 
of acp (AS = 0) and acp (AS = 1) are circles in the (p, fj) -plane. These constraints are 
complementary to the information from other observables in K and B physics. 

4) Inclusive direct CP-asymmetries are theoretically clean: The uncertainties can be con- 
trolled and systematically reduced by higher order calculations. 

5) The CP-conserving observables Br (AS = 0), Br (AS = 1) and Br NL (B — > no charm) 
define circles in the (p, rf) -plane centered on the p-axis. 

6) Br (AS = 0) is well suited to determine \ V u b/Vcb\, with little sensitivity to 7. 

7) "Double penguin" contributions, which are part of the next-to-next-to-leading order, en- 
hance Br(AS= l)by 13%. 
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6 TEN MESSAGES FROM THIS WORK 




Figure 5: The black annulus shows the region allowed by a measurement of Br (AS = 0). The 
dark shading corresponds to Btnl(B — > no charm) and the lightly shaded area shows the con- 
straint on (p, rj) obtained from Br (AS = 1). 



8) The present incomplete next-to-leading order (NLO) result imposes a large p-dependence 
on Br (AS = 1). Here a calculation of all NLO corrections to penguin operator matrix 
elements is necessary. 



9) Br (AS = 1) exceeds Br (AS = 0) by roughly a factor of two. 



10) The determination of \V u b/V c b\ from Br NL (B — ► no charm) is competitive to the standard 
method from semileptonic decays, once the NLO calculation mentioned in 8) has been 
done. 
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Figure 6: The lightly shaded area shows the constraint stemming from acp (AS = 0) and the 
dark shading marks the area allowed from acp (AS" = 1). 
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A. Exact formulae 

Here we show how the quantities entering sect. 3 and sect. 4 are related to the results of [10]. 
These expression are useful for readers who are not satisfied with the approximate formulae 
for F, L B , L a , etc. They are also helpful, if one wants to calculate the branching ratios and 
CP-asymmetries in extensions of the Standard Model. Then one needs to change the Wilson 
coefficients entering the T^ 's accordingly. 
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A EXACT FORMULAE 



Now F [14] reads 

F = 



G 2 F ml\V cb \ 2 
64vr 3 r tot 



1 + 



2ml 



64tt 3 



1 + 



SL 



2m? 



h{xl) [1 + a s ^)/{2n)h SL {x c )] - 6 (1 - xlf X 2 /m 2 b 



(29) 



in terms of the notation of [10]. Here we have used the common trick to evaluate T tot = T S l/Bsl 
via the semileptonic rate and the experimental value of the semileptonic branching ratio B SL . 
This eliminates various uncertainties associated with the theoretical prediction for T tot . The non- 
perturbative corrections involving the kinetic energy parameter Ai has been factored out in (29), 
because A x cancels in Br and A C p. 

Likewise for the decay rates corresponding to the quark level transition b — > qqd', q = u,d,s and 
d! = d, s, one has 



r = 

-■- nil. 



t £ cfii 

i,j=l 




a. 



r*u — 



tc 



Ttt = 




+ 5b ij + -^2Re [hij + gniO)] 

47T 



+ t5b ij + ^(g ij (0)+tg* i (x c )) 



a. 



C$C 2 — b 2 %t 
lix 



+ 5k 



+ ^-C s £ C s b sa 

27F i=3,..6 



47T 



C\ k 22 ^x c , x c , 



A 

m b 



(30) 



Here t = 1 for b — > uud', while t = for b — > ssd' and 6 — > ddcf . The C/s, a s and the loop 
functions hij, gij in (30) are understood to be evaluated at the scale \x = 0(m b ). The T^'s depend 
sizeably on /i and x c as indicated in the approximate formulae in sect. 3. Further they depend 
on m t and M w , this dependence, however, is marginally small. g 22 (x c ) = g(x c , n/m b ) is the 
fundamental penguin functions entering all g^'s. For this work we have newly calculated 



#42 X, 



JL 



#62 Xc, 



A 

m b 



932 x, — 
m b 



m b 



g 0,-^ )+g 1,-^ 
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m b 

27 



n f g 0, 



/i 



A 



+ 5 U c , 



m b 



952 



A 



+ 9 1 



A 



a:, — = 0, 
m 6 , 



...... 98 8 32 2 

"^ ln - + y-7! 7r+ ^ 7r - 



71/ = 3, 



(31) 



m b 8 y/3 

These quantities correspond to the diagrams of Fig. 4 with q' — u, q — u, d, s, c, b and the left cut 
marking the final state uud. For the remaining g^'s we refer to [10], where also analytic formulae 
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for g(x,fi/m b ) and the ft^'s and hi/s [17] can be found. In (30) the leading nonperturbative 
corrections are also included, the Sbij's [18] depend on A 2 = 0.12 GeV 2 . The values in (31) 
correspond to the NDR scheme, the vanishing of g 52 involves in addition the standard finite 
renormalization of Q§ introduced in [19] and related to the definition of the "effective" coefficient 

C 8 . 

Another new result is T cc in (30). We have calculated the "double penguin" contribution stem- 
ming from the square of Fig. 2 with q' = c. Although being of order a 2 this term is numerically 
relevant in AS* = 1 decays, because it is proportional to Cf and the tree-level result is CKM 
suppressed. We have also included T cc in the AS = coefficients of (17). Approximately one 
finds 

k 2 2 Oc, x c , Xp) = (l - [l.52 - 11.5 (x c - 0.3) + 7 (x c - 0.3) 2 

+ (l.84- 5.6 (x c - 0.3) - 19(x c -0.3) 2 )ln^ + 0.791n 2 xJ.(32) 

with r = for the quark final states uud, uus, ssd, dds, and r = 1 for ddd and sss. The result 
in (32) receives corrections of order (x c — 0.3) 3 and reproduces k 2 2 with an error of 2.6 % for 
0.25 <x c < 0.35and 0.5 < fi/m b ^2.0. 

Our results for Br (AS = 0) and Br (AS = 1) also include the decay rates for b — > s g and 
b — > dg. Here to order a s all Tj/s are zero except for 

T tt (b^sg) = Y tt (b^dg) = ^^C 8 2 . 

The approximate formulae in (17), (21) and (24) further correspond to a s (M z ) = 0.118 cor- 
responding to a s (fi = 4.8GeV) = 0.216. The dependence on a s (M z ) is non-negligible, but 
smaller than the //-dependence. 

When calculating A C p for the inclusive AS = or AS = 1 final state, we must add the r^'s 
for b — > uud', b — > ssd' and b — > ddd'. Imr (u contains lmg(0), which, however, cancels 
when summing Imr tu for the three decay modes b — > uud', b — > ssd' and b — > ddd', so that 
A C p(AS = 0) and A C p(AS = 1) vanish for x c > 1/2 as required by the CPT theorem. The 
cancellation takes place when summing the contributions of different cuts of the diagrams in 
Fig. 4 as found in [2]. 

One comment is in order here: The terms of order a s in (30) depend on the renormalization 
scheme. This originates from the fact that when renormalizing H in (1) one already uses the 
unitarity relation £ u + £ c + £ t = 0. After using this relation to eliminate, say, £ c in (2) one 
finds the coefficients of |£ u | 2 , |£ t | 2 and £ t £* scheme independent. Consequently by changing the 
scheme one can shift terms in A C p in (7) from e.g. the term proportional to sin 7 to the one 
multiplying sin (3. This scheme ambiguity, however, is suppressed by a factor of C^q/Ci^ with 
respect to the dominant contribution to a C p. The constraints on (p, rj) derived from Br and a C p 
are scheme independent, of course. 
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